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We compute the graviton one loop contribution to a classical energy in a traversable wormhole 
background. The form of the shape function considered is obtained by the equation of state p = up. 
We investigate the size of the wormhole as a function of the parameter u. The investigation is 
evaluated by means of a variational approach with Gaussian trial wave functionals. A zeta function 
regularization is involved to handle with divergences. A renormalization procedure is introduced and 
the finite one loop energy is considered as a self-consistent source for the traversable wormhole. The 
case of the phantom region is briefly discussed. 



I. INTRODUCTION 

The discovery that our universe is undergoing an accelerated expansion [1| leads to reexamine the Friedmann- 
Robertson- Walker equation 

a 47T 

- = -^-(p + 3p), (1) 
a 6 

to explain why the scale factor obeys a > 0. Indeed, it is evident from the previous formula, that a sort of dark energy 
is needed to cause a negative pressure with equation of state 

P = up. (2) 

A value of lu < — 1/3 is required for the accelerated expansion, while w = — 1 corresponds to a cosmological constant. 
A specific form of dark energy, denoted phantom energy has also been proposed with the property of having lu < — 1 . 
It is interesting to note that the phantom energy violates the null energy condition, p + p < 0, necessary ingredient 
to sustain the traversability of wormholes. A wormhole can be represented by two asymptotically flat regions joined 
by a bridge. To exist, it must satisfy the Einstein field equations: one example is represented by the Schwarzschild 
solution. One of the prerogatives of a wormhole is its ability to connect two distant points in space-time. In this 
amazing perspective, it is immediate to recognize the possibility of traveling crossing wormholes as a short-cut in 
space and time. Unfortunately, although there is no direct evidence, a Schwarzschild wormhole does not possess 
this property. It is for this reason that in a pioneering work Morris and ThorneQ and subsequently Morris, Thorne 
and YurtseverJU studied a class of wormholes termed "traversable" . Unfortunately, the traversability is accompanied 
by unavoidable violations of null energy conditions, namely, the matter threading the wormhole's throat has to be 
"exotic" . It is clear that the existence of dark and phantom energy supports the class of exotic matter. In this direction, 
LoboQ, Kuhfittig[5[ and SushkovQ have considered the possibility of sustaining the wormhole traversability with 
the help of phantom energy. In a previous work, we explored the possibility that a wormhole can be sustained by 
its own quantum fluctuations 7]. In practice, it is the graviton propagating on the wormhole background that plays 
the role of the "exotic" matter. This has not to appear as a surprise, because the computation involved, namely 
the one loop contribution of the graviton to the total energy, is quite similar to compute the Casimir energy on a 
fixed background. It is known that, for different physical systems, Casimir energy is negative and this is exactly one 
of the features that the exotic matter should possess. In particular, we conjectured that quantum fluctuations can 
support the traversability as effective source of the semiclassical Einstein's equations. However in Ref.0], we limited 
the analysis in the region where the equation of state® assumes the particular value lu = 1. In this paper, we 
will consider lu S (0, +00) [H, although the semiclassical approach can be judged suspicious because of the suspected 
validity of semiclassical methods 1 at the Planck scale @. The rest of the paper is structured as follows, in section HT1 
we define the effective Einstein equations, in section iLUl we introduce the traversable wormhole metric, in section ITVl 
we give some of the basic rules to perform the functional integration and we define the Hamiltonian approximated up 
to second order, in section fVl we study the spectrum of the spin-two operator acting on transverse traceless tensors, 
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1 To this purpose, see also paper of Hochberg, Popov and Sushkov llOU and the paper of Khusnutdinov and Sushkov[llfl. 
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in section PVTl we regularize and renormalize the one loop energy contribution and we speculate about self-consistency 
of the result. We summarize and conclude in section IVTlI 

II. THE EFFECTIVE EINSTEIN EQUATIONS 

We begin with a look at the classical Einstein equations 

Gfi V — kT^ v , (3) 

where is the stress-energy tensor, G Mt , is the Einstein tensor and k = 8irG. Consider a separation of the metric 
into a background part, g^, and a perturbation, h^ Vl 

g^=g^ + h^. (4) 

The Einstein tensor can also be divided into a part describing the curvature due to the background geometry 
and that due to the perturbation, 

G^v {gap) = G^y {(jail) + AG^ (<? Q /3, h a [j) , (5) 

where, in principle AG M „ (g a /3) ^a/j) is a perturbation series in terms of h^. In the context of semiclassical gravity, 
Eq.® becomes 

G M „ = k (T^y en , (6) 

where (T^„) e " is the renormalized expectation value of the stress-energy tensor operator of the quantized field. If the 
matter field source is absent, nothing prevents us from defining an effective stress-energy tensor for the fluctuations 
as 2 

(T^) ren = (AG M „ (g a p, h a0 )) ren . (7) 

K 

From this point of view, the equation governing quantum fluctuations behaves as a backreaction equation. If we fix 
our attention to the energy component of the Einstein field equations, we need to introduce a time-like unit vector 
such that u ■ u — — 1. Then the semi-classical Einstein's equations ([6]) projected on the constant time hypersurface £ 
become 

G„„ {g aP ) uV = k (T^u») ren = - (AG M „ (g aP , h a p) u»u") ren . (8) 
To further proceed, it is convenient to consider the associated tensor density and integrate over E. This leads to 3 

/ d 3 x^G^ (g aP ) u»u v = - f d 3 xH^ = --L / Sx^g (AG^ (g a07 h a0 ) u»u v ) ren , (9) 

where 

H (o) = i^kl _ J_^(3) (10 ) 

V 9 

is the background field super-hamiltonian and Gijki is the DeWitt super metric. Thus the fluctuations in the Einstein 
tensor are, in this context, the fluctuations of the hamiltonian. To compute the expectation value of the perturbed 
Einstein tensor in the transverse-traceless sector, we use a variational procedure with gaussian wave functionals. In 
practice, the right hand side of Eq.Q will be obtained by expanding 

4 0) + 4 1} + 4 2) + 



_<*lff E |*>_\* 

^wormhole ~ 



(ii) 



2 Note that our approach is very close to the gravitational geon considered by Anderson and Brill ll3| . The relevant difference is in the 
averaging procedure. 

3 Details on sign conventions and decomposition of the Einstein tensor can be found in ApeendixlBl 
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and retaining only quantum fluctuations contributing to the effective stress energy tensor. H$ represents the hamil- 
tonian approximated to the i th order in and $ is a trial wave functional of the gaussian form. Then Eq.@ 
becomes 



4 0) = / d 3 xH^ = 



The chosen background to compute the quantity contained in Eq.© will be that of a traversable wormhole. 



(12) 



III. EINSTEIN FIELD EQUATIONS AND THE TRAVERSABLE WORMHOLE METRIC 

In Schwarzschild coordinates, the traversable wormhole metric can be cast into the form 

dr 2 



dsf = -exp(-2<£(r))dr 



2 XjcU i „:„2 a .2 



[<W 



sin 



(13) 



i 6(r) 
r 

where (j> (r) is called the redshift function, while b (r) is called the shape function. Using the Einstein field equation 

G>„ = kT^, (14) 
in an orthonormal reference frame, we obtain the following set of equations 

1 V 



p(r) 



8irG r 2 ' 



(15) 



Pr (r) 



8irG 



-ii 

r \ r 



b(r) 



(16) 



n 1A b(r) 



in i it i jj 



b'r-b 



(17) 



in which p (r) is the energy density, p r (r) is the radial pressure, and pt (r) is the lateral pressure. Using the conservation 
of the stress-energy tensor, in the same orthonormal reference frame, we get 



Pr = ~ (Pt - Pr) ~(P + Pr) 4>' ■ 



The Einstein equations can be rearranged to give 



b' = 8ttG P (r) r 2 , 



(18) 



(19) 



, _ b + 8TrGp r r 3 

Now, we introduce the equation of state p r = up, and using Eq. fTB"]) , Eq. P0|) becomes 

b + ub'r 



2r 2 ( i _ Mr) 



The redshift function can be set to a constant with respect to the radial distance, if 

b + ujb'r = 0. 

The integration of this simple equation leads to 



(20) 



(21) 



(22) 



(23) 
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where we have used the condition b(r t ) = r t . In this situation, the line element (fl"3|) becomes 

dr 2 

ds 2 = -Adt 2 + + r 2 [dO 2 + sin 2 8d V 2 } , (24) 

l-(^) + " 

where A is a constant coming from <j)' — which can be set to one without loss of generality. The parameter to is 
restricted by the following conditions 

6 > t )<i; ^-o. =► (25) 



T r— »+oc 



Proper radial distance is related to the shape function by 



where the plus (minus) sign is related to the upper (lower) part of the wormhole or universe and where 2F1 (a, b; c; x) 
is a hypergeometric function. Two coordinate patches are required, each one covering the range [rt, +00). Each patch 
covers one universe, and the two patches join at rt, the throat of the wormhole defined by 

r t = min{r(0}. (27) 
When u> = 1, we recover the special case where b (r) = r 2 /r, then the line element simply becomes 

ds 2 = -dt 2 + + r 2 [ d6 2 + sin 2 ed ^ (28) 



and Eq. pS|) simplifies into 



dr' 



l(r) = ± ^= = ±Jr 2 -r 2 => r 2 =l 2 +r 2 t . (29) 
^ a/1-4 



The new coordinate I covers the range — 00 < I < +00. The constant time hypersurface E is an Einstein-Rosen bridge 
with wormhole topology S 2 x R . The Einstein- Rosen bridge defines a bifurcation surface dividing E in two parts 
denoted by E + and E_. To concretely compute Eq . fHi]) . we consider on the slice E deviations from the wormhole 
metric of the type 

9ij = gi 3 + hij, (30) 
where gij is extracted from the line element (|24|) whose form becomes 

ds 2 = -dt 2 + /,,,</.,■' </r< . (31) 

IV. ENERGY DENSITY CALCULATION IN SCHRODINGER REPRESENTATION 

In order to compute the quantity 

r d?x^m {AG^ (g af3 , h af >) u»u v ) ren , (32) 

we consider the right hand side of Eq. (|12p . Since is linear in hij and h, the corresponding gaussian integral 

disappears and since 

^G^igativPu" = -H, ( :-!:-!) 
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it is clear that the hamiltonian expansion in Eq. IT2|) does not coincide with the averaged expanded Einstein tensor of 
Eq. ({3"2")) because Eq. (|33[) involves a tensor density. Therefore, the correct setting is 



d 3 xy^g 



(34) 



where ^ is the second order expanded tensor density weight. Following the same procedure of Refs.[l5l[l6l]. the 

potential part of the right hand side of Eq. ((54"]) becomes 



/ d 3 x^g --hAh + h H Ah H - \h^ViVih\ + l-hViV.h 11 - \wR ia h a 3 + \hRi_ 
J 4 4 Z Z Z Z 



//'■- + -m,' ! ii, 



-hRh 



The term 



/ d 3 xy/g 



-Rh i: >h vj - -hRh 
4 J 4 



(35) 



(36) 



makes the difference between the hamiltonian expansion and the Einstein tensor expansion. To explicitly make 
calculations, we need an orthogonal decomposition for both -k^ and hij to disentangle gauge modes from physical 
deformations. We define the inner product 



(h, k) := / y/gGpHhy (x) k kl (x) d 3 x, (37) 

by means of the inverse WDW metric Gyiy , to have a metric on the space of deformations, i.e. a quadratic form on 
the tangent space at h, with 



QijH = (gikgjl + g il g jk _ g^ij g kl ^ 

The inverse metric is defined on co-tangent space and it assumes the form 



(38) 



(p, q) := / s //C /,,.,//< (x) q kl (x) d 3 X, 



so that 



(39) 



(40) 



Note that in this scheme the "inverse metric" is actually the WDW metric defined on phase space. The desired 
decomposition on the tangent space of 3-metric deformations [13, EH is: 



hij = g%j + ( L 0ij + h ij 



(41) 



where the operator L maps £j into symmetric tracefree tensors 



(iO <J - = v i c- + v^i--5y(V-o. 

Thus the inner product between three-geometries becomes 

(h, h) := f x gC'^'h,, (x) h kl (x) d 3 x 



(42) 



V9 



(43) 
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With the orthogonal decomposition in hand we can define the trial wave functional 

9 {hij 0?)} = M exp {--L UhK-'h)^ + ((L£) K- 1 (L£))l „ + (hK^h) 



Trace 
x,y 



(44) 



where J\f is a normalization factor. We are interested in perturbations of the physical degrees of freedom. Thus we 
fix our attention only to the TT tensor sector reducing therefore the previous form into 



* {hij (1?)} = AAcxp {-- (hK-'h^l 



(45) 



Therefore to calculate the energy density, we need to know the action of some basic operators on 9 [hij]. The action 
of the operator hij on |\&) = 9 [hij] is realized by 



h tj (x) |*) = h tj (~x)iS>{hij} . 



(46) 



The action of the operator 7Tjj on \9), in general, is 



TTij (x) |*) 



Shij fx ) 



The inner product is defined by the functional integration: 

(*i I *2> = J [Dhij] 9\ {h tJ }9 2 {h kl } , 



(47) 



(48) 



and by applying previous functional integration rules, we obtain the expression of the one-loop-like Hamiltonian form 
for TT (traceless and transverse) deformations 



1 



ijkl (16vrG) 



(A 2 r^ {x,x) l 



ikl 



The propagator K 1 - {x^x) iakl comes from a functional integration and it can be represented as 



K 1 - (I?.-?) 



takl 



E 



2A(r) 



(49) 



(50) 



where h^ X (~x) are the eigenf unctions of A 2 , whose eigenvalues will be denoted with E 2 (r). t denotes a complete 
set of indices and A (r) are a set of variational parameters to be determined by the minimization of Eq. (|49]) . The 
expectation value of H is easily obtained by inserting the form of the propagator into Eq. 



(16ttG) A, (t) + 



(16ttG) A, (t) 



(51) 



By minimizing with respect to the variational function A^ (r) , we obtain the total one loop energy for TT tensors 



E 



TT 



E 2 (r) + jE 2 (r) 



(52) 



The above expression makes sense only for Ef (r) > 0, i — 1,2. The meaning of Ef will be clarified in the next 
section. Coming back to Eq. (fl"2"|) , we observe that the value of the wormhole energy on the chosen background is 4 



f d 3 xH (0) = _ 1 f d 3 xVdR (3) =A[u) n 



lj > — 1 



(53) 



4 Details of the calculation can be found in the Appendix fX] 
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where 

A H = — - — B (\ -±-) = //M. (54) 

B (x, y) is the Beta function and T (x) is the gamma function. Then the one loop the self-consistent equation for TT 
tensors becomes 

A^) r -±=-E TT . (55) 



Note that for the special value of u> = 1, we get 



= -£ TT > (56) 



in agreement with the result of Ref. [7j. Note also that the self-consistency on the hamiltonian as a reversed sign 
with respect to the energy component of the Einstein field equations. This means that an eventual stable point for 
the hamiltonian is an unstable point for the effective energy momentum tensor and vice versa. 



V. THE TRANSVERSE TRACELESS (TT) SPIN 2 OPERATOR FOR THE TRAVERSABLE 
WORMHOLE AND THE W.K.B. APPROXIMATION 



In this section, we evaluate the one loop energy expressed by Eq. (|52p . To this purpose, we begin with the operator 
describing gravitons propagating on the background Q24p . The Lichnerowicz operator in this particular metric is 
defined by 

{^h TT )\ := - {& T h TT ){ + 2 (Rh TT )\ - R {h TT )\ , (57) 
where the transverse-traceless (TT) tensor for the quantum fluctuation is obtained by the following decomposition 

ti = hi-lti*+l%h=(h T )\ + ±fi>h. (58) 

This implies that (/i T )* = 0. The transversality condition is applied on (/i T )^ and becomes V 'j (^ T )^ = 0. Thus 

(A T h^)i = -As(h^y i + ^(l- h -^), (59) 



where A 5 is the scalar curved Laplacian, whose form is 



b(r)\ d 2 /4r~b'{r)r~3b(r)\ d L 2 



r 



A -=( 1 -^)^+ ^ -J*-^ (60) 



and is the mixed Ricci tensor whose components are: 

Ra _ f b'(r) b(r) y(r) b(r) V (r) 6(r)| 
1 \ r 2 r 3 ' 2r 2 2r 3 ' 2r 2 2r 3 J ' [ ' 

The scalar curvature is 

b' (r) 

R=m= 2 -^ (62) 

We are therefore led to study the following eigenvalue equation 

(A 2 h TT ) j .=E 2 hj (63) 

where to 2 is the eigenvalue of the corresponding equation. In doing so, we follow Regge and Wheeler in analyzing 
the equation as modes of definite frequency, angular momentum and parity [l9j. In particular, our choice for the 
three-dimensional gravitational perturbation is represented by its even-parity form 
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h% m (r,0,<l>)=diag 



H(r)[l 



b(r) 



, r 2 K (r) , r 2 sm z §L (r) 



with 



H{r) = h\ (r)-lh(r) 
K(r)=h 2 (r)-\h(r) 
L(r) = hUr)-lh(r) 



(64) 



(65) 



From the transversality condition we obtain h\ (r) = h\ (r) . Then K (r) = L (r) . For a generic value of the angular 
momentum L, representation (f6"4"]) joined to Eq. (T5"9")) lead to the following system of PDE's 



-A, + 2 - b P ~ H (r) = E 2 U H (r 



(- A, + 2 + ^ - ^)) i^r) = ^Jf (r) 



where A; is 



, _|. | b(r)^ d 2 + (4r-b'(r)r-3b(r)^ d 1(1 + 1) 6 ^ b(r) 



Defining reduced fields 



dr 2 \ 



H(r) 



2r 2 



AM. 



dr 



K{r) 



Mr) 



r r 

and passing to the proper geodesic distance from the throat of the bridge defined by 

dr 



dx = ± 



6(r) 



the system (f66|) becomes (r = r (x)) 



(66) 



(67) 



(68) 



(69) 



where we have defined r = r (x) and 



with 



-£ + V 1 (r)\f 1 (x)^ElJ 1 (x) 
-£ + V 2 (r)}f 2 (x)=ElJ 2 (x) 

Vi (r) = + Ui (r) 

V 2 (r) = M + f/ 2 ( r ) ' 



C/i (r) - Cl (r) + (i - 3) c 2 (r) , 
C/ 2 (r)= Cl (r) +3(1 + 1) c 2 (r) , 



ci (r) 



(70) 



(71) 



(72) 



In order to use the WKB approximation, we define two r-dependent radial wave numbers ki (x, I, Ei jTl i\ and 
k 2 (x,l,E 2 ^ 

k 2 (x,l,^)=% M -mp.- Ul {T) 



k\ {x,l,E 2 . n i 



(73) 



Elm ~ ^ - U 2 (r) 
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The number of modes with frequency less than Ei, i = 1, 2, is given approximately by 

g(Ek) = J Vi(l,Ei) (2l + l)dl, (74) 

where Vi (l,u>i), i — 1,2 is the number of nodes in the mode with (l,uii), such that 

1 f + °° I 

Vi(l,Ui) = — dxJkf (x,l,LJi). (75) 

27r J-oc V 

Here it is understood that the integration with respect to x and I is taken over those values which satisfy kf ^x, I, Eij > 
0,i = l, 2. Thus the total one loop energy for TT tensors is given by (recall that r = r (x)) 



where 



E 



TT 



1 2 r+oo dg i Ei) 2 r+oo I" i r+oo , 



/■+C» 

1 ' 



drr 2 [pi + p 2 ] , 



Pi 



Pi 



4tt X 



j_ r+ 00 

4t Jy/U 2 (x) 



E*JEl-U 2 (r)dE 2 



(76) 



(77) 



VI. ONE LOOP ENERGY REGULARIZATION AND RENORMALIZATION 

In this section, we proceed to evaluate the one loop energy. The method is equivalent to the scattering phase shift 
method and to the same method used to compute the entropy in the brick wall model. We use the zeta function 
regularization method to compute the energy densities p\ and pi. Note that this procedure is completely equivalent to 
the subtraction procedure of the Casimir energy computation where zero point energy (ZPE) in different backgrounds 
with the same asymptotic properties is involved. To this purpose, we introduce the additional mass parameter p in 
order to restore the correct dimension for the regularized quantities. Such an arbitrary mass scale emerges unavoidably 
in any regularization schemes. Then we have 



Pi (e) 



1 

47T 



2e 



+ 3o 



dEU 



Ef 



(E? U { (r)) 



(78) 



If one of the functions Ui (r) is negative, then the integration has to be meant in the range where Ef + Ui (r) > 0. In 
both cases the result of the integration is 5 



64tt 2 



hi 



P 



Ui(r) 



21n2- 



1 



(79) 



where the absolute value has been inserted to take account of the possible change of sign. Then the total regularized 
one loop energy is 



E TT (r t ,e;n)=47t {2 I dr- 

>rt 4 /i_Mri 



[(Pi (e)+P2 (£))] 



(80) 



where the factor 4-7T comes from the angular integration, while the factor 2 in front of the integral appears because 
we have come back to the original radial coordinate r: this means that we have to double the computation because 
of the upper and lower universe. To further proceed, it is useful to define the following coefficients: 



5 Details of the calculation can be found in the Appendix |Pl 
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1. 



and 



2. 



a := 2r t I dr 



1 - 



6(r) 



[U?(r) + U$(r)] 



b := - 



+00 

2r t I dr 



1 - 



6(r) 



[^(r)ln|^(r)| + C^(r)ln(U- 3 (r))] 



(81) 



(82) 



c := 2r t / dr 



b(r) 



[Ul (r) + XJ\ (r)] In (p 2 ) = 2alnp. 



Rescale the radial coordinate r so that p — r/rt, then 



6(r) /r t \i+i /i 



(7) 



r \ r / \p 
Thus the coefficient a simply becomes 

/+00 
dp 

while the coefficient b changes into 



; tfi (r) = E/x (p) /rf; *7 2 (r) = C/ 2 (p) jr\. 



? [U?(p) + Ui(p)], 



1 - 



(») 



1 



(83) 



(84) 



(85) 



+00 2 r+00 2 

6 := -2 / dp g [C/f (p) In \U t (p)| + E/f (p) In (C/ 2 (p))] +2 / rip g [17? (p) + [/f (p)] (2 lnr t ) 



1 



/+00 
dp 



1 - 



1+ 



= [J7 1 2 (p)ln|i7 1 (p)| + C/ 2 (p)ln(C/ 2 (p))] + 2alnr t - b + 2a\nr t . 



The result of the integration over the r coordinate leads to the following expression 



a b 2a I \^8rtfi 
sr t r t r t \ \fe 



and the self consistent equation ()55p can be written in the form 



*<»>5-i5F 



a b 2a I V^rtU \ 

1 1 ln ±r 

sr t r t r t \ ye J 



(86) 



(87) 



Actually, also the coefficients a and b depend on lo, but at this level is not relevant. In order to deal with finite 
quantities, we renormalize the divergent energy by absorbing the singularity in the classical quantity. In particular, 
we re-define the bare classical constant G 



1 



G G eA (w) 167rr t 2 
Therefore, the remaining finite value for the effective equation ([57]) reads 

rt 1 



A(u) 



G 16tt 



n n 



(88) 



(89) 
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This quantity depends on the arbitrary mass scale It is appropriate to use the renormalization group equation to 
eliminate such a dependence. To this aim, we impose that[2l| 



d (A{uj)r t 



dfi \ G (//) 



d_ I _L_ 

^ dfj, I 16n 



namely 



A {to) nil 



dGp 1 ( M ) 
9/i 



87rr t 



= 0. 



(90) 



(91) 



Solving it we find that the renormalized constant Go should be treated as a running one in the sense that it varies 
provided that the scale \i is changing 



1 1 K\ ( [l 



or G (p) = 



G (Mo) 



K 



l + GoM^ln(^)' " ^MSTrrf 
where no is the normalization point. We substitute Eq. (f9"2")) into Eq. (j59")l to find 



1 



4j 

G {no) 167T 



b 2a ( V8r t no 
12 + 72 ln -^7T- 



where we have divided by r t . In order to have only one solution 6 , we find the extremum of the r.h.s. of Eq. 
we get 



a - b _ ( V8r t no 
2a ~ n I #e 



and 



anl 



G (no) 2-kA (u) y/e 
With the help of Eqs.(JSU» and |95|) . Ea.((9"2] ) becomes 



exp 



a — b 



exp 



a — b 
2a 



In 



GoO) Go (mo) A{uj)8irr? \noJ A {uj) 2-K^fe 



/' 



a/4 



exp 



a — b 



1 + 2 In 



Mo 



(92) 



(93) 
and 

(94) 
(95) 



1 



1 + 2 In | ± 

Mo 



G (Mo) 

It is straightforward to see that we have a constraint on n/fj-o -Indeed we have to choose 

fj, > no exp ( — —) = .6065306597/ito; 



(96) 



(97) 



otherwise Go (/i) becomes negative . We have now two possibilities: 1) we identify Go (^o) with the squared Planck 
length, then the wormhole radius becomes 



n 



' 4j 

WirA (w) 



(98) 



Note that in the paper of Khusnutdinov and Sushkov lll| . to find only one solution, the minimum of the ground state of the quantized 
scalar field has been set equal to the classical energy. In our case, we have no external fields on a given background. This means that 
it is not possible to find a minimum of the one loop gravitons, in analogy with Ref. Moreover the renormalization procedure in 
Ref. is completely independent by the classical term, while in our case it is not. Indeed, thanks to the self-consistent equation J55}, 
we can renormalize the divergent term. 
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where we have reestablished the u> dependence of the coefficient a. It is useful to write the expression for to — > +00 
and for u> — ► 0. We get 



f(l + Mli-2 1n(2))+OK 2 )) 



— > +00 



(99) 



^7= + O (^2) 



L w — > 



The following plot shows the behavior of r t as a function of u>. It is visible the presence of a minimum for u> — 3.35204, 




FIG. 1: Plot of the wormhole throat r t as a function of u in the positive range. 



where f t (lu) — 1.11891. As we can see, from the expression (|99p and from the Figfl] the radius is divergent when 
oj — > 0. At this stage, we cannot establish if this is a physical result or a failure of the scheme. When oj — > ±00, 
ft approaches the value 1.15624/ p , while for ui = 1, we obtained ft = 1. 15882Z p . It is interesting to note that when 
lu — > +00, the shape function b{r) in Eq. (|23j) approaches the Schwarzschild value, when we identify f t with 2MG. In 
this sense, it seems that also the Schwarzschild wormhole is traversable. 
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2) We identify /j,q with the Planck scale and we get from Eq. (|94|) the following plot Note the absence of a minimum. 




CO 



12 3 



Differently from case of Eq. (|98p . there is no access to the phantom range. In our case, we suppose that the graviton 
quantum fluctuations play the role of the exotic matter and even if we fix the renormalization point at the Planck 
scale as in Ref.|ll|. we find a radius f t > r w . 

VII. SUMMARY AND CONCLUSIONS 

In this paper, we have generalized the analysis of self-sustaining wormholesQ by looking how the equation of state 
(|2"|) can affect the traversability, when the sign of the parameter u> is positive. The paper has been motivated by the 
work of Lobo[l|, Kuhfittig[f| and SushkovQ , where the authors search for classical traversable wormholes supported 
by phantom energy. Since the phantom energy must satisfy the equation of state, but in the range us < —1, we 
have investigated the possibility of studying the whole range — oo < us < +oo. Unfortunately, evaluating the classical 
term we have discovered that such a term is well defined in the range — 1 < us < +oo. The interval — 1 < us < 
should be interesting for the existence of a "dark" energy support. Once again, the "dark" energy domain lies outside 
the asymptotically flatness property. So, unless one is interested in wormholes that are not asymptotically flat, i.e. 
asymptotically de Sitter or asymptotically Anti-de Sitter, we have to reject also this possibility. Therefore, the final 
stage of computation has been restricted only to positive values of the parameter us. In this context, it is interesting 
to note that also the Schwarzschild wormhole is traversable, even if in the limiting procedure of us — ► +oo. Despite of 
this, the obtained "traversability" has to be regarded as in "principle" rather than in "practice" because the wormhole 
radius has a Planckian size. We do not know , at this stage of the calculation, if a different approach for a self sustained 
wormhole can give better results. On the other hand, the positive us sector seems to corroborate the Casimir process 
of the quantum fluctuations supporting the opening of the wormhole. Even in this region, we do not know what 
happens approaching directly the point us = 0, because it seems that this approach is ill defined. Nevertheless, in this 
paper we have studied the behavior of the energy. Work in progress seems to show that dealing with energy density 
one can get more general results even in the "phantom" sector [8[. 

Acknowledgments 

The author would like to thank S. Capozziello, F. Lobo and the Referee for useful comments and suggestions. In 
particular, the author would like to thank the Referee, for having brought to his attention the paper of Ref.[6|. 



14 



APPENDIX A: COMPUTATION OF THE CLASSICAL TERM 



Here, we give details leading to Eqs. (|53l 54|) . We begin with the definition of the hamiltonian which in the static 
case simplifies into 



J d 3 xH^ = J d 3 xV§R {3) - (Al) 



With the help of Eqs (|23l 24|) , we get 



'f*p) s ' I:=1 + I. (A2) 

Gio 1^ \ r ) If \s a; 



In the previous integral there is an extra factor "2" coming from the counting of the universes. We change variable 
to obtain 

2r t f + °° dt _ r t ( ' 1 1 \ _ r t 

G^Jo cosh 2 " 2 /- (t) ~G(T+Tj) U' " G A H (A3) 

with 











/5F r 




(1- 











In Eq. (|A3p , we have used the following formula 



cosh" t 2 V 2 ' 2 7 \ Re (/x - i/) < v; 



APPENDIX B: EINSTEIN EQUATIONS AND THE HAMILTONIAN 

Let us consider the Einstein equations 

G^v = R^v - \^9p.vR = k8ttGT^ v . (Bl) 

R^v is the Ricci tensor and R is the scalar curvature. If it M is a time-like unit vector such that g^u^u" = — l,then 
the Einstein tensor G^ becomes 

G^u" = R^v? - X -g^ l u v R = R^ l u v + Ir. (B2) 

By means of the Gauss-Codazzi equations[22], 

R = i? (3) ± 2R pu u>*u v tK 2 ± K^K^, (B3) 

where is the extrinsic curvature and R^ is the three dimensional scalar curvature. For a time-like vector, we 
take the lower sign and Eq. (|B2p becomes 

G^uW = \ (i? (3) + K 2 - K^K^) . (B4) 
If the conjugate momentum is defined by 

^0, = V — y , Rgl j, v _ K ^^ ( B5 ) 
2k 
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then 



K 2 - K^K^ 



2k 

7W, 



(B6) 



and 



/Wq 



2 k 



2k ( 7T 



/(3)g V 2 



(B7) 



namely Eci. tflT))) with the reversed sign. 



APPENDIX C: ENERGY COEFFICIENTS 

In this appendix, we report details on computation of the coefficients of the one loop graviton energy. We begin 
with the coefficient a defined by 



a := r t / dr 



a/1 



b(r) 



[Ul (r) + Ul (r)] , 



(CI) 



which, with the help of Eqs. ([72]) , becomes 



r+oo 

a = r t I dr 

Jr t 



2c{ (r) + 2 ( 9 + - + 4 (r) + - Cl (r) c 2 (r) 



(C2) 



where 



6 /, 



Mr) = 3 1- 7 



(C3) 



Defining the dimensionless variable p = we can write 



oi = ?*< / 



+ OC 



, ^ ■ j N (r)] = 72 r dp U^)-l) % p ->-^) 



l + w V 2 1 + oj 



(C4) 



a,2 = n I dr- 

* ^1-(^) J 



■ ^> • - • 4 ) • <'•) 



-if. " " 



!^)r d "(" (,+i) - i )" i ^" §(,+i) 



1 / 6 _5_\ cj / 1 3cj + 2 

2 1 + cj + cj 2 /1 + cj U'l + w 



(C5) 



and 



03 = n I dr 



1- a 



-ci (r) c 2 (r) 



24 /-+ 00 , / 



( 1 +^)_l N )%-2-f(l+i) 



24 /3 2w + l 



1 + uj V 2 l + w 



(C6) 
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so that 



a = a (u) — 



1+OJ 



72B 



5 u> 



2' 1+uW 2 



1 /„ 6 5 \ „ / 1 3u; + 2\ 24 „ fi 2uj + 1 
9+ - + — I B I -,- I + — B 



2' 1+w / u \2' 1 + w 



where B (x, y) is the beta function.. The same procedure applies to the coefficient b. 



/+oo 
dp 



(i) 



1+ 



= [Uf(p)ln\U 1 (p)\ + Ul(p)ln(U 2 (p))]. 



We can separate the logarithmic functions 

ln|f7i(p)| =ln 12w(p( 1+ ^) - l) + (1 - 3w) -ln(2w)- (s+^J 111 / 3 



6i c 



ln|J7 2 (p)| =ln 12w(p( 1+ i) -l)+3(l+w) -ln(2w)- (3+^Jlnp 



b2a 



b-2b 



&2c 



and compute separately the different terms. fe 16 + 6 2 b simply becomes 

&ii + fob 



2m(2u;) J dp — t = In (2a;) a (a; J . 



-(i) 



The coefficient &i c + 62c is 



foe + foc= 3 + 



=)r* 



p 2 [^(p) + ^ 2 (p)] 







iA-( 





ln/9, 



namely 



3 + 



1 



(j) 



2c? (p) + 2 ( 9 + - + A ) c 2 ( p ) + - Cl ( p ) C2 ( p ) 



w or 



In p 



(C7) 



(C8) 



(C9) 



(CIO) 



(Cll) 



(C12) 



The first term, b/,i will be summed with the coefficient a, while the second term can be written as bf^i + &/,62 + &/,&3 
with 



2 ,+QO 



^ 00 d / ;(p( 1 +^)-l)%- 2 -i( 1 +^)lnp 



72 3 + 



1 \ / co \ I ,3 



Tu + 5 

dss 5 (1 + s) 2 < 1+ "» ln(l + s) 



1 



72 3 + - 



w/ ll+w/ \2 1 + oj 



2(1+ a;) 



1 + w 



(C13) 



W/ V W U J Jl I /,\l+77 



1 



(i) 
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o 1 \ 1 fn 6 5 

3+ _ 9+ _+ 

OJ I A \ OJ 0J Z 



( 3+ ^)K 9+ ^)(x^) / 



2 r + OO 

OJ \ f 1 7 " + 5 

' ' dss~i(l + s) 2 <!+-) ln(l + s) 



and 



2(l + w) 



1 + w 



3 + 



OJ / OJ 



-0) 



1+7 



OJ OJ 



^ +0 °dp(p( 1+ =)-l)%- 2 -i(^)lnp 



(C14) 



/ 1\ 24 / W \ 2 /" + 00 1 7^ + 5 

(3 + -j-( TT ^jy o ^(l + -r^m(l + -) 



= 3 + 



1 \ 24 / w 



W/ W \ 1 + w 
where we have used the following relation 

-iln (l + x) 



B 



3 2oj + 1 

2' 1+OJ 



7w + 5 
2(l + w) 



•(t£) 



r+oo 

Jo ' 



dxx^ 



/ = (M, !/-/*)[* M - * (i/ - /i) + In 7 ] , 



(7 + x) 

where ^ is the digamma function. Globally, the coefficient b\ c + b 2c becomes 



(C15) 



(C16) 



= 3 + - 



oj I V 1 + w 



725 1* " 



2' 1 + w 



2(1+0;)' 



1 + w 



24 /3 2w + l 
H i> 



w \2' 1+w 



Ioj + 5 
2(l + w) 



2cj + 1 



1 + OJ 
(C17) 

The term &i a + b 2 b is more complicated. The logarithmic term is of the form In (7^2 + Sx). It is useful to re-express 
it in terms of hypergeometric function, namely 



with 5 = \2oj and 



In (71,2 + Sx) = S— 2 Fi [1,1; 2; -5— ) + In 712, 

71,2 V 71,2/ 



71 — 1 — 3oj 
72 = 3 (1 + oj) 



(C18) 



(C19) 



Expression 



/+00 
dp 



1 - 



0) 



1+77 



becomes 



In 7! 



OJ 



1+OJ 



365 



5 oj 



2' 1 + oj 



'■!(/>)+ ("J- 3 ) ^(P) + 2(^-3)n (/))ro(» 



oj J \ 2' 1 + oj J \oj J \ 2' 1+w 



(C20) 



(C21) 
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and expression 



/+oo 
dp 



0) 



1+7 



becomes 



In 72 



1+LJ 



36S(i " 



2 1 + cj 



<j(p) - (^ + 3j ^(^ + 2^+3)0 (,)],,(,;] 



3 3w + 2\ /3 „\„/3 2w + l 



(C22) 



(C23) 



Concerning the hypergeometric part, after having defined the variable s = , we use the following formula 

f 

Jo 



dx- — : — -p 2F1 (a, b; c; — ujx) dx = z a ~ p B (a, p — a) 3F2 (a, b, a; c, a — p + 1; uiz) 



la (x + zf 



_ a T(c)T(a-a + p)T(b-a + p)T(a- p) 
T(a)T(b)T(c-a + p) 



3 F 2 (a - a + p, b - a + p,p; c- a + p,p - a + 1; uiz) , (C24) 



which, in our specific case, becomes 



~ — J dsj^—^pS 2F1 ( 1, 1; 2; — <5— — ] ds = B (a, p - a) 3 F 2 ( 1, 1, a; 2, a - p + 1; - — 



7i,2 



71,2 



6 y a T{l-a + p)T(l-a + p)T(a-p) 



,71,2/ T{2-a + p) 

Therefore for 71 we get 

S /■+«, 2^(1,1; 2; -l(pH)_l)) p 2 
7i 7i 



3 F 2 1 - a + p, 1 - a + p, p; 2 - a + p, 1 - a + p; . (C25) 

71,2, 







I— 1 





- 1) 



<A (p) + ( - - 3 J C2 (p) + 2 I - - 3 ) ci (p) c 2 (p) 



Every piece leads to: a) 

^ r dP , p2cl{p) 1 - 1) a * r 1, 1; 2; -— (p(-^ - 1) 

71 7i / 77\i+I v ' V 71 v ' 



(C26) 



= 36— (^Y^J J dss * (! + s) _,TTT ^ 2*1 (^1, 1; 2; -— s 



7i 



ft) 



+ 



— (- 

71 V w 



36 A 



71 \ 1 + Ul 



1 2' 1 + W 



r ('"•-•) , ' - 
— 3-^2 



3 F 2 [1,1,-; 2,— — ; — 
1 w — 1 7lu + 5 2u w — 1 5 



l + w'l + w'2(l + w)'l+w'l + w'7i 



5 
71 



(p(^i) - 1) 



(C27) 
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6li 



+00 



4w (1 + uj) J 



1 $ 
dss2(l + s) 2 Fi ( 1, 1; 2; s 



<$7i 





2u + 1\ 




HI- 


1 + OJ ) 









2« + l 


r 


f 3w+2 X 


; 


2 / 
r ( 


2u+l\ 
1+c J 


r 


/ 4w+3 
v 1+w 


) 



3-F2 



71 / 4w (1 + u) 
/3w + 2 3w + 2 7w + 5 4w + 3 3w + 2 <S 

\l + £j' 1 + W '2(1+0;)' 1+Cj' l+Cj'71 



and 



1 + lo 71 _ 



(C28) 



6(5 r°° 3 7^ + 5 / ,5 

/ dss* {l + s y^+^ 2 Fi 1,1; 2; s 

1+Wo V 71 



_65_ 

1 + CJ 



2' 1 +uj 



■ ;F, ! i.l,|;2, ' 



2' ' 1+w' 71 



7i, 



r 


f 2w+l> 




2 






r ( 


3w+2 S 
1+w 


) 



3-F2 



/2w + l 2w + l 7cj + 5 3w + 2 2w + l <5 
V 1 + w ' 1 + w ' 2(1 +w) ! 1 + uj ' 1 + uj ' t7 



(C29) 



The same set of expressions works for 72 too. 



APPENDIX D: THE ZETA FUNCTION REGULARIZATION 

In this appendix, we report details on computation leading to expression (|78[) . We begin with the following integral 

/+ = M 2£ Jo + °° duj- 



P(e) 



(u 2 + U(x)f-2 



(uJ 2 -U(x))'-2 



with U (x) > 0. 



(Dl) 



1. 7+ computation 



If we define t = u>/ y / ^U~(x), the integral I + in Eq. (|D1[) becomes 



P {e)=^U{xf~ 



+00 .2 
dt 



(i 2 + 1) £ 



1 -^U{xf-* B (\,s-2 



W(^ r(i)r(£ - 2) -^— 2 



r( £ -i) " 4 ^Vw; r( £ -|)' 

where we have used the following identities involving the beta function 

r+00 t 2x ~ 1 

B(x,y)=2 dt x- 

V ,J> Jo (t^ + l) x+y 

related to the gamma function by means of 

T(x)T(y) 



M 2 V r(e-2) 



Rex > 0, Rey > 



B(x,y) 



T(x + y) 



(D2) 



(D3) 



(D4) 
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Taking into account the following relations for the T-function 

r(i + e ) 



r( £ -2) = 



e(e-l)(s-2)- 



r e- 



i\ r( B + i) 



£ — 



and the expansion for small e 



rn f.-i = i- :.: + (>•(.--)■ r( e + i] =rQ] - e rQ) ( 7 + 2in2) + o(e^ 



(D5) 



x £ = 1 + elnx + O (e 2 ) , 



where 7 is the Euler's constant, we find 



p{e) = - 



U 2 (x) 
16 



e 



r/(V; +21n2 4 



(D6) 



(D7) 



2. 7_ computation 



If we define t = u>/ \JU (x), the integral 7_ in Eq. (|Dl[) becomes 



p( £ )= M 2£ (7(x) 



\2-e 



+00 



- ' l =l^U(xf^B(e-2, 3 --s 

{t 2 - ip 3 2 V 2 



/ ^( x) >-.E3- e ) r < e - 2 > 1 — 2 ( ^ 



4^F 



u( x y 



where we have used the following identity involving the beta function 



Ib(i-v-H,v 
V \ P 



+oc 



dtt"- 1 (t p -lf 



p > 0, Re v > 0, Re fi < p — p Re v 



and the reflection formula 

r (z) r (1 - z) = -zv (-z) r (z) 

From the first of Eqs. (|D5P and from the expansion for small e 

r (|- £ ) = r (T) (l- £ (-7- 2 m2 + 2)) + 0( £ 2 



(D8) 



(D9) 



(D10) 



x £ = 1 + elnx + O (e 2 ) , 



we find 



U 2 (x) 
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1 ' 1:1 v,< 



2 In 2 - - 
2 



(Dll) 



(D12) 
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